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Abstract 

A family of local models containing two angles as hidden variables 
is denned for experiments measuring polarization correlation of optical 
photons. Searching for the best model of the family, that is giving 
predictions most close to quantum mechanics, allows deriving Bell- 
type inequalities which may be tested with relatively low detection 
efficiency. 

PACS. 03.65.Ud Entanglement and quantum nonlocality (e.g. 
EPR paradox, Bell's inequalities, GHZ states, etc.) - 42.50.Ar Photon 
statistics and coherence theory 

1 Introduction 



More than forty years have elapsed since John Bell[l] proposed his celebrated 
inequalities. These inequalities, which involve measurable quantities, provide 
necessary conditions for local realism and, in some experiments with ideal 
set-ups, contradict the predictions of quantum mechanics. Many empirical 
tests have been performed of local realism against quantum mechanics, via 
the Bell inequalities, but no experiment has been conclusive. In fact although 
the results have generally agreed with the predictions of quantum mechan- 
ics, no experiment has given results incompatible with local realism, as is 
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shown by the existence of local realistic models for all of them[2]. The in- 
ability to perform a true empirical test of local realism is commonly disguised 
with the claim that it has already been refuted by the experiments, modulo 
some irrelevant loopholes. But I think that the extreme difficulty to make a 
loophole- free test, proved by the unsuccessful effort of forty years, does not 
support the common wisdom that the question of local realism is settled. On 
the contrary the conclusion is that further research is needed. 

Locality should be understood in the relativistic sense, that is as im- 
possibility of superluminal communication. Thus testing locality requieres 
measurements made in regions spatially separated, in the sense of relativ- 
ity theory. As a consequence the tests are extremely difficult with massive 
particles [3] and reliable experiments must be performed with photons. On 
the other hand, tests with high energy photons are not possible due to the 
lack of efficient polarization analyzers. These difficulties have caused that 
most of the experimental tests have been performed with optical photons, 
where good polarization analyzers exist and locality may be insured. How- 
ever these experiments suffer from a detection loophole, due to the fact that 
a good overall detection efficiency has not yet been achieved. In fact, it is 
well known that efficiencies as high as 80% are required for loophole-free tests 
of local realism. Detectors with high quantum efficiency already exist, but 
there are other difficulties reducing the overall efficiency to about 30% or less 
in practice. 

It is common to formulate the question of local realism in terms of local 
hidden- variables (LHV) theories. That is, an experiment refutes local realism 
if there is no LHV model compatible with their results. Bell inequalities 
are necessary conditions for LHV theories but, as said above, they are very 
difficult to test. Actually all inequalities empirically violated till now are not 
genuine Bell inequalities, derived from the conditions of realism and locality 
alone, but inequalities whose derivation requires auxiliary assumptions. It 
is rather obvious that violations of such inequalities do not refute the whole 
family of LHV theories but only restricted families, namely those fulfilling the 
auxiliary assumptions. In my opinion the families of local hidden variables 
theories so far refuted by the experiments are rather unplausible[2] . This 
is the case, in particular, for those fulfilling the fashionable "fair sampling 
hypothesis" . This is the assumption that the photons actually detected are 
representative of the whole set of photons emitted. But the role of hidden 
variables is precisely to distinguish, from each other, several physical systems 
in the same pure quantum state, that is systems which are identical according 
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to the standard interpretation of quantum mechanics. Thus in any LHV 
model of an experiment it is natural to assume that the photons detected 
and those not detected correspond to different values of the hidden variables 
and consequenly the sample of detected photons is not representative of the 
whole set. In conclusion, the fair sampling assumption amounts at dismissing 
all sensible hidden variables theories from the start. 

In the present paper I study a new family of LHV models which is rather 
natural, in my opinion, and it allows the derivation of inequalities able to dis- 
criminate between the said family and quantum mechanics. The inequalities 
are easily testable in optical experiments, in particular the tests would require 
only moderate detection efficiencies, of the order of 30%. The said family 
was already considered in a previous publication [4], where an inequality ful- 
filled by some sub-family was proposed. That inequality has been tested 
empirically with the result that it was fulfilled and the quantum predictions 
contradicted although, according to the authors of the experiment [5], the 
contradiction cannot be considered a violation of quantum mechanics. In 
the present paper I propose inequalities which should hold true for all mem- 
bers of the said LHV family, but are contradicted by the quantum predictions 
in some cases. 

2 A natural family of local hidden variables 
models 

For the sake of clarity I shall consider experiments measuring the polarization 
correlation of optical photon pairs, although the generalization to other cases 
is possible. The set-up consists of a source of photon pairs each member of 
the pair travelling in a different direction, crossing a lens system, a polar- 
ization analyzer and arriving at a detector. If the polarization planes of the 
analyzers are determined by the angles and 2 , respectively, the results of 
the experiment may be summarized in two single rates, Ri(4>i ) and i?2(02 ) ? 
and a coincidence rate Ru((pi , 02 )• ( in recent experiments four coincidence 
rates, rather than one, are measured because two-channel polarizers are used, 
a situation which will be considered below). In a polarization correlation ex- 
periment the detection rates should be obtained from appropriately defined 
probabilities pi,p2 and pu, that is 

R A<f>j ) = R 0Pj(<t>j ), 3 — 1? 2, #12(01 , 02 ) = RoPni^i , 02 ), (1) 
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where Rq is the production rate of photon pairs in the source, a quantity not 
measurable in standard experiments. (Here it is assumed that all photocounts 
come from photons produced in pairs in the source. In practice there may be 
counts of a different origin, e.g. dark counts in each detector, but they will be 
neglected for the moment, see section 4 below). Following Bell a LHV model 
consists of three functions, /(A), Qi(X, Qz(X, 4>2)i where A stands for one 
or several hidden variables, such that the detection probabilities could be 
obtained by means of the integrals 



PM,) = J f{\)QiMi)d\ Pi2(0!,0 2 ) = J /(A)g 1 (A,0 1 )g 2 (A,0 2 )dA. 

(2) 

The essential requirements of realism and locality imply that the said func- 
tions fulfil the conditions 

/(A) > 0, J f(X)d\ = 1, < Qj(\, ^) < 1. (3) 

A natural, but relatively simple, family of local hidden variables model 
is obtained if we assume that the set A of hidden variables may be written 
A = {xi, fa, Xii P2} j where Xi an d Pi (X2 an d A* 2 ) are variables of the first 
(second) photon of a pair and Xj is a polarization angle, so that Xj an d 
Xj + 7T represent the same polarization. Actually we may assume that fi l and 
/i 2 label a set of variables each, rather than a single one. In principle the 
function /(.) gives the correlation amongst the four (sets of) variables. But 
simplifies a lot the model to assume that the variables, or sets of variables, 

and fi 2 are uncorrelated amongst themselves (this is a consequence of 
locality) and uncorrelated with the polarization angles (this is an assumption 
of simplicity), so that the functions /(A) and Qj(X, -) may be written 

/(A) = p(xi, x 2 )gi M 92 M , Qj(\ 4>j) = QAxj,Pj, <Pj), (4) 

where p and gj are positive and normalized. When eqs.(4) are inserted in (2) 
and the integrals in ■ performed we get 

Pj{<t>j) = J p(Xi,X2) p j(XjAj)dXidX2, 

Pl 2 (01,0 2 ) = J P(Xl,X2) P l(Xl,0l)P2(X2i ( f ) 2)dXldX2, (5) 
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where 

Pj(Xj, <f>j) = J 9j(»j)Qj(Xj, /';' fyWj- ( 6 ) 

Thus our model requires defining only the functions p and Pj, which ful- 
fil conditions of positivity and normalization similar to those of / and Qj, 
respectively, in eqs.(3) . 

In order to derive testable inequalities from our LHV model we shall 
consider firstly experiments where there is rotational symmetry, which is 
common at least approximately in actual experiments. That is, we assume 
that the single rates Rj do not depend on the orientation of the polarizers 
(i. e. the angle -) and the coincidence rate depends only on the difference 
0i — 02 = 0- We shall also assume symmetry amongst the polarizer-detector 
systems so that P 1 = P 2 . (See section 4 below for the study of cases where 
these symmetries do not hold true). In these conditions it is appropriate to 
replace eqs.(5) by 

PM) = J p(Xi - X2) p {Xi - <t>i)P{X2 - <t>2)dXidX2i (7) 
Pj = J p(Xi - X2)P{Xj - (f>j)dXidX2i 3 = 1, 2. (8) 

where here and below all functions are periodic with period tt and the in- 
tegrals go from -7r/2 to tv/2. In addition the functions p and P possess the 
following properties of positivity, symmetry and normalization (p is normal- 
ized so that p 12 {4>) = 1 if P = 1 in (7)) 

pM^(-,)>o,/^^i/,,o<p W ^(-,)<i. (9) 

I shall add the following two conditions which are plausible on physical 
grounds 

d 4^<o, ^<o. (10) 

d\x\ d\x\ 

The first inequality means that the pairs where \\ = X2 are m ost likely 
produced in the source (this equalitity may correspond to parallel or per- 
pendicular polarization, depending on the actual experiment). The second 
inequality means that the detection is most probable when the incoming pho- 
ton has the polarization close to the plane of the analyzer. Eqs.(7) to (10) 
define our "natural" family of LHV models. 
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For later convenience I make the change of variables 

Xi - 0i = u, x 2 ~ 02 = v > 0i - 02 = 0, (11) 

which leads to 

p 12 (0) = J p( u - v + cj))P(u)P{v)dudv. (12) 
Hence we get the following results 

Pi 2 (0)# = P1P2, Pi = P2 = Co/tt, y Pl2(0) cos 20 rf0 = Cf/V, (13) 
where the constants are defined by 

C k = J P{x) cos (2kx) dx. (14) 
It is also convenient to introduce the new function 

f(y) = f(-y) = Jp(x + \)P{x - y -)dx. (is) 

whence we get 

Pi 2 (0) = J P (y + 4>)f{y)dy = J P (y)f(y - 4>)dy, (16) 

the latter equality following from the periodicity of the functions involved. 

3 Comparison with the quantum predictions 

Now I will investigate whether the quantum predictions are compatible with 
the family of LHV models above defined. Quantum mechanics predicts, for 
experiments with rotational symmetry and both detectors having the same 
efficiency, 77, 

P? = \v, P12 (0) = \rf (1 + V cos 20) . (17) 

Our family of models (7) agrees with the quantum-mechanical prediction for 
the single probabilities, pj, provided we choose 

Co = J P(x)dx = ^. (18) 
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Then I will search for the best LHV model of the form (7), defining "best" 
by the condition that, for fixed rj and V, the prediction for the coincidence 
probability is as close as possible to the quantum prediction. That is the 
quantity S must be a minimum, where 



S = 



Pi 2 (0)-^ 2 (1 + Vcos20) 



(19) 



with pi2((f>) given by (16) and (15) and the functions p(x) and P(x) fulfilling 
the conditions (9) . For the solution of the problem I shall proceed in two 
steps. Firstly, for a given P(x) fulfilling (18) , I search for the best positive 
and normalized function p(x). In the second step I shall look for the best 
P(x). 

In the first step fixing P(x) amounts at fixing the function f(y) (see 
eq.(15)). For the solution of the variational problem, eq.(19) , the positivity 
of p is insured writing 

p(x) = ip(x) 2 , ^eR, 

and the normalization constraint may be taken into account by means of a 
Lagrange parameter A . Thus the problem becomes 



5 | S [^(x) 2 ] + \j tjj(x) 2 dx | = 0. 

Hence we get, after some algebra, that either ip(x) = or 

j ^{xfdx J d<Pf(x-^f(y-^ = \ + ^r ] 2 J #/(y-0)[l + ycos20]. (20) 

For the solution of this integral equation the properties of p(x), see (7) , 
suggest the Fourier expansion 



^{xf = J2 Ak cos ( 2kx )i ( 21 ) 
which, inserted in eq.(20) leads to (see eq.(14)) 

MC\ cos(2ky) = X + -if [C 2 + VCl cos(2y)] . (22) 



k=0 
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Here I have used the following relation, easily derivable from the properties 
of the function f(x) (see eq.(15)), 



f(u — v) cos(2ku)du = C\ cos(2kv). 



Eq. (22) holds true for any y if 
4A + rfC 2 



A, 



4C 4 



A 1 



rfV 
4C?' 



for k > 2. 



Thus the function p (x) will be 

4A + rfC 2 



p(x) 



4C% 



1 + 



v 2 c*v 



C 2 (4A + ^C 2 ) 



cos 2x 



(23) 



(24) 



+ 



where [.] + means putting if the quantity inside the parenthesis is negative 
and A should be calculated from the normalization condition (see eq.(9)). 
There are two cases which must be analyzed separately. 
The first case corresponds to the function P(x) fulfilling 



C\ f P(x) cos 2x dx i— 

— = — — > vV 

C JP(x)dx 



Then we may take A = leading to 



7T Z 



c 2 

1 + -±V cos 2x 

C\ 



7l z 



ir 2 ri 2 
1 + —^V cos 2x 
AC{ 



(25) 



(26) 



where the choice of Co, eq.(18), has been taken into account. This gives 
perfect agreement with the quantum predictions, (17) , for both single and 
coincidence probabilities. 

It is interesting to know the range of values of i] and V where it is possible 
the agreement between the model and quantum predictions. Actually the 
constraints (9) and (18) put an upper bound to the left hand side of the 
inequality (25) . It is not difficult to realize that the best function P(x), that 
is the one saturating the bound, is 



This choice leads to 



P(x) = 1 if \x\ < ixrj/A, otherwise. 

d= 5111^/2), 



(27) 
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whence (25) becomes 

V < V m 



sin 2 (nrj/2) 
H/2) 2 



ir 2 r] 2 



(28) 



the latter approximation being valid for rj << 1. We stress that, if this 
inequality is fulfilled, there are LHV models of the type (7) giving complete 
agreement with quantum mechanics. In fact, the choice eq.(26) for p (x) and 
eq.(27) for P(x) gives the desired result. 

We pass to analyze the second of the two cases in eq. (24) which, by the 
arguments of the previous paragraph, will correspond to the violation of 
the inequality (28) . In this case the function p(z) may be written, correctly 
normalized according to eq.(9) , 



1 r / % i-i / cos2z\ / 

p (z = - 7r + tan (2e ) - 2s] 1 + — ,£6 0, 

7r y cos 2e J , V 



7T 



(29) 



and I will search for the departures between the model and quantum-mechanical 
predictions using that function. For the single detection probability, pj, there 
is no departure, that is the predictions exactly agree, provided we make the 
choice (18). However there is disagreement for the coincidence probability, 
which may be most conveniently exhibited expanding (29) in a Fourier series 
of the form 



(30) 



n=0 



We get a — 1, as it should p (z) being normalized in the sense of (9), and 



ai 



[tt + tan(2e) - 2e\ 
2 



7T 



2e 



+ sin (2s) 



n(n z 



1) 
16 



;-iy 



cos (2e 

sin (2ne) — n tan (2e) cos (2ne) 
tt + tan(2e) - 2e 



l + 2e 2 



7T 



+ 0(e 4 ), 
(31) 



= (-l) n — e 3 + 0(e 5 ) forn> 2, 
3n 

Hence it is straightforward to get the model predictions for p i2 from (12) and 
(27), that is 



P12 



n=l 



r<2 



cos(2n<; 



(32) 
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where C n was defined in eq.(14) . 

Our aim is the find the best local model, in the sense of S (19) being 
a minimum for given 77 with P(x) and p(z) fulfilling, respectively, eqs.(18) 
and (29) . It is not difficult to realize that the minimum S will correspond to 
Pi2(0) being of the form 

P12 (0) = ^ 2 [1 + V cos (20) + 5 (0)], (33) 

with 5 (<f)) containing only terms with n > 2 in the Fourier expansion (32). 
For any choice of P(x) this will be the case if the following relation between 
e, 77, V and C\ holds true 

y= -S*r- (34) 

Now we shall choose P[x) so that |5(0)| is as small as possible. From 
eqs.(31) we see that this requires that e is small which, as a\ increases with 
e, implies that C\ must be high (see eq.(34)). Hence the best choice for 
P(x) will correspond to the maximum possible value of C\ compatible with 
eq.(18) ,which happens if P(x) is chosen as in eq.(27). After that the value 
of e may be obtained from eq.(34) giving, to order 0(e 2 ), 



an equality clearly showing that e is a measure of the violation of the in- 
equality (28) (e = if the inequality holds true). 

With the said choices of P(x) and p (x) the departure of the model from 
the quantum predictions is given by 

5(0) = £a n ^^cos(2n0). (36) 

Inserting here the expressions for a n , eq.(31) , it is straightforward, although 
lengthy, to get the expression of 5 (0) . However the parameter e will be very 
small in practice (i. e. e « tt/4), which allows working to lowest order in e, 
that is 0(e 3 ) (see (31)). Although that approximation of a n is good only for 
the terms with small n, the terms with high n contribute but slightly. We 
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may get the explicit form of 5 (0) to order 0(e 3 ) from eq.(36) with a n given 
by eq.(31). We shall use the summation formula 



E 

n=l 



:-ir 



X 2 TT 2 



cos (nx) — — , X G [— 7T, 7r] , 



12' 



but we must substitute (— 1)™ cos(n(x — n)) for cos(nx) whenever x > it (but 
x < 27r) and a similar change if x < —n. In this case the summation formula 
needed is 



n=l 



I \ 2 2 

. . (X — 7T) 7T . „ . 

cos (nx) = - - — , x e [0, 2tt] 



Thus the term S ((f)) is given up to order 0(e 3 ) by 



^in 2 {rcrj/2) 
1 (^/2) 2 



cos (20) 



2 

1 + - 

7] z 



v + -101-1 



7T 



©(l0|-f(l-^)) 



(37) 

where 6 (x) is the Heavside function, fulfilling (x) — 1 if x > 0, © (x) = 
if x < 0. (Of course we shall take <5 (0) =0 if eq.(28) holds true.) 



4 Testable inequalities for the family of local 
models 

A parameter giving a quantitative measure of the discrepancy between quan- 
tum mechanics and the best model of the family defined by eqs.(7) and (8) 
is, from eq.(36) , 

71 J 1 (nnri/2) 



The summation in n would be straightforward but lenghty and, as all terms 
in the sum are positive, a lower bound is provided by the first few terms. 
With just one term we get 

A > y/2sin 3 (2£) sin 2 ^) ^ 8^2 sin 2 (nrj) =3 

~ 3 [(tt — 2e) cos2e + sin2e] 7r 2 i] 2 37r n 2 rf 
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Actually it is easy to calculate A to order 0(e 3 ) either performing the sum 
involved in eq.(38) or directly integrating the square of the expression (37). 
I get 




1 
2 



sin 4 (nrj/2) 
H/2) 4 



V - 



sin 



(^/2)\ 3/2 _ 



h/2) 2 ; 



= D(rj) : (40) 



where I have taken e from eq.(35) (as before (.) , means putting zero if 
the quantity inside the bracket is negative). The quantity D{rj) gives the 
deviation, as defined in eq.(37) , between the best local model and quantum 
mechanics. Thus in any experiment whose results are compatible with the 
family of local models defined by eqs.(7) and (8) , the deviation will be greater 
than D(rj). That deviation A cxp may be obtained from the measurement of 
the coincidence detection rates, i?i 2 (0), at different angles, taking 8 (0) to 
be the difference of that quantity appropriately normalized minus the best 
cosinus fit of Pi2(0)- Thus we get 



*cxp 



V cos 20 



1/2 



> D(ri), 



(41) 



whilst quantum mechanics predicts A cxp = . This inequality has the virtue 
that it may be tested even at relatively low detection efficiencies (provided 
that V is close enough to unity as is usual), so removing the biggest obstacle in 
the standard tests of Bell's inequalities. It replaces the previously proposed 
inequality (14) of Ref.[4], which is stronger but valid only for a family of 
LHV theories more restricted than (7). 

From the practical point of view the inequality (41) requires the mea- 
surement of the coincidence detection rate at different values of the angle, 0, 
between the polarizers. I may consider n angles defined by 



3 



n 



71, j = 1,2.. .n, 



and approximate the integral in the left hand side of (41) by a sum, so leading 
to the inequality 



■'cxp 



n 

n ^ 



R 



121 

JR 



1 - V cos 20, 



12/ 



1/2 



> D( V ), (42) 
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where 

3=1 

The quantity V should be taken from the quantum prediction for the 
actual experiment, but the inequality holds true for any value of V provided 
we use that same V in the calculation of D(r)), eq.(41) . It is interesting to 
use the value of V giving the minimum value to A exp for the given empirical 
rates {Ru((j)j)} , that is 

_ o E? = i^i2(0 J )cos20 J 
where I have taken into account the equalities 

n n 

^cos20 J =O,^cos 2 20 J = ^. 

3=1 3=1 

In this case I get 



A min 




(44) 

this inequality making possible empirical tests of the family of local models 
defined by eqs.(7) and (8) without any reference to quantum mechanics, 
whilst eq.(42) allows tests of the local models versus quantum mechanics. 

The test of the inequalities involve the measurement of n/2 + 1 rates if n 
is even or (n + l)/2 rates if n is odd, once we take into account the symmetry 
Ri2((f>j) = #12 (?r — 4>j). For even n > 4, eq.(44) may be related to the Fourier 
expansion eq.(38) using the interpolation formula 

^1 = 1 + ^6,008 (2^.). ( 4 5) 
whence the left hand side of (44) may be interpreted as 

1/2 

I 1 '"/^ \ 



1 min 



n/2-1 ) ^ 

\Y. h l + h l/A • ( 46 ) 

k=2 J 
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The most simple non-trivial case corresponds to n — 4, that is the choice 
of the angles 0, 7r/4 and ir/2. In this case eqs.(41) and (44) lead to the 
inequality 

A . = *■»(«) + Ms) -"Mi) > D( ) (47) 

Rn (0) + Rn (f ) + 2Rn (f) ~ 1 1 ' 

Strictly speaking, from eq.(44) we may only conclude that the modulus of 
the left hand side fulfils the inequaltiy, but we may take that side as positive 
because 02 > in the expansion (30). At a difference with (41) the inequality 
(47) is not rigorous due to the approximation made in going from eq.(41) to 
(42) , the possible error being greater as the number of measured rates is 
smaller. 

In addition to the inequalities it is possible to derive some predictions of 
the proposed local models for measurable quantities. For instance, defining 

_ gl2 (0) - gig (f ) _ r- Rl2 (f) -Rl2 (f ) , . 

Va = MoKMf)' Vb = V2 i? 12 (f) + i? 12 (f)' (48) 

their difference may be estimated. The former quantity is called the visibil- 
ity (or contrast) of the polarization correlation curve and the latter is the 
quantity measured in the standard tests of Bell's inequalities. According 
to quantum mechanics both quantities should be equal to the parameter V, 
(17), but in our model there is a difference. In fact, from eq.(45) it follows 
that 

= ; 1 r , r - 61 + 63 - 6 2 - 64, Vb = — ~ 61 - 6 3 + 64, 

I + 62 + 64 1 - 64 

where I have taken into account that, from a comparison of the empirical 
quantities bk with the model quantities of eqs.(36) and (35) , we get 

64 <N <6 2 «1~6 1 = V. 
Hence the family of LHV models here studied predicts (see eqs.(31) and (35)) 

V B - V A K h - 2 h + 2bi ^ K (v- V2 , (49) 

3tt V. (irri/2) J + 

K being a number of order unity. In some experiments it is reported, instead 
of Va, the value of V, which is obtained from a fit to the empirical coincidence 



14 



rates R 12 ((f>j). In this case the model predicts 

I stress that these estimates are correct only for the "best" local model of the 
family. Thus their value is only indicative of the order of magnitude to be 
expected in actual experiments if they are compatible with the said family 
of local models. 

In recent optical tests of Bell's inequalities people use two-channel po- 
larizers and four coincidence rates are measured, p ++ (0) , p + _ ((f)), p_ + (0) , 
p (0) . Typically there is symmetry between the two channels, at least ap- 
proximate. In this case our model may be extended to these experiments 
using two functions P + ((f)) and P_ (0) = P + ((f) + tt/2) instead of only one, 
P ((f>) , as in eq.(7) . Assuming the form (27) for P + ((f)) we get the prediction 

((f) + tt/2) = R- + ((f) + tt/2) = R ++ ((f)) = R- ((f)) , (51) 

the latter given by eq.(33) (putting e — if the inequality (28) holds true). 
The quantity reported in the experiments is the correlation, defined by 

F U] = (0) + R- (0) - R+- (0) - R-+ (0) , v 

W) R ++ ((f)) + R__ ((f)) + R + _ (0) + R_ + ((f)) ' 1 ] 

for which our model predicts 

E((f) — V cos(20) — 5 (tt/2 + 0) , (53) 

where 5 (0) is given by eq. (37) . The parameter most accurately measured in 
the experiments is 

S =\3E (tt/8) - E (3tt/8) | = 2y/2V B , (54) 

which is predicted to have the value 2v^2 according to quantum mechanics 
in ideal experiments with 100% detection efficiency, but S should be less 
than 2 for any LHV model fulfilling the "fair sampling hypothesis" (see the 
Introduction section.) 

From eq.(48) it is easy to see that, in our LHV family of models, the 
quantities Vb and Va defined in (48) correspond, respectively, to 2v^2 times 
the quantity S of eq.(52) and to 

V A = \ (E(0) - E(n/2)) . (55) 
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This allows testing the model prediction eq.(49) in experiments involving 
two-channel polarizers from measurements of the correlation E(cf>), eq.(52) , 
at the four angles 0, 7r/8,37r/8 and tt/2. Also the prediction eq.(50) may be 
tested getting Vb from (54) and V from a fit of E(0) as shown in eq.(53) . 



5 From ideal to real experiments 

In the previous sections I have analyzed rather idealized experiments. In 
actual polarization correlation experiments involving parametric down con- 
verted photons the results obtained are not so simple as assumed in eqs.(l) 
and (17) . Indeed it is frequent that the coincidence detection rate is not ro- 
tationally invariant in the sense that R\ 2 (</>!, <t> 2 ) does depends on the angles 
4> 1 and 02 separately rather than on the difference = <p l — 2 only. In 
addition, the single rates may depend on the positions of the polarizers, that 
is Rj = Rj((f)j), and the efficiencies, r] 1 and r] 2 , of the detectors may not be 
the same. Finally not all photons may arrive in pairs at the detectors (ether 
because there is some single photon production or because a fraction of the 
photons do not enter the apertures.) The relevance of these nonidealities will 
be analyzed in the present section. 

The quantities which may be measured in typical experiments involving 
single channel polarizers are: Two single rates, R\ and R 2l the rate R\ 2 
of counts in the second detector conditional to counts in the first detector 
(within an appropriate time window) and the quantum efficiencies of the 
two detectors, i] 1 ,i] 2 (these measured in auxiliary experiments). In order 
to apply the analysis of Bell[l] to the experiments we should define single 
and coincidence probabilities. But for the test of the family of LHV theories 
studied in this paper (defined by eqs.(7) and (8)) we do not need the single 
probabilies and the single rates, Rj, are not used. In order to obtain the 
coincidence probability, p\ 2 , I introduce the rate, Ro, of pair production in 
the source, whence I may get the p\ 2 using Bayes rule of probability theory 
as follows 

Ri R12 R12 / r „N 
Pn = Pi x p 2/1 = — x — = — . (56) 
rto rC\ Ko 

The adequacy of this definition of p\ 2 requires that all photons leaving the 
source are produced in pairs, with no more than one pair within each time- 
window. That is, eq.(56) amounts at neglecting single photon production, 
dark counts and the possible production of several photon pairs within one 



16 



time-window. I shall assume that all these difficulties are avoided by a sub- 
traction of "accidental coincidences" , which is a standard practice in polariza- 
tion correlation experiments. (I stress that this is so because we are attemting 
to test a restricted family of LHV models, given by eqs.(5) , but would not 
be valid if we attempted to refute the whole family of LHV theories [2].) Ac- 
tually the pair production rate i? cannot be measured in the experiments 
but this is not a real difficulty. Indeed the inequalities fulfilled by the LHV 
family of models proposed in this paper involve only the dependence of the 
coincidence probability on the positions of the polarizers, determined by the 
angles <p 1 and 2 . In the experiments this dependence is usually studied by 
measurements with a polarizer at a fixed position, say 2 fixed, an the other 
polarizer's position varied. Typically the measurements are made at two 
different positions of the second polarizer, say <p 2 — and 2 = 7r/8. An 
appropriate interpolation leads to the following parametrization 

pi 2 a R 12 oc 1 + [V m + (V M - V m ) cos40 2 ] cos 20, (57) 

where Vm(Kw) is the maximum (minimum) visibility of the detection coinci- 
dence curve when the angle + is varied. 

A local hidden- variables model able to reproduce eq.(57) for not too high 
detection efficiencies is given by the functions p and Pj (see eq.(5)) 

p{Xi,X 2 ) = ^ {1 + [W m + (W M -W m ) cos Ax 2 ] cos 2 X }, 
Po{x v <t> 3 ) = J = 1,2, (58) 

with any /5 - e [rj j: l] and appropriate choices for W M and W rn < W M - 
Here \ = Xi ~ X2 an d @ (•) * s the Heavside function (x) — 1 if x > 0, zero 
otherwise. However if the detection efficiencies are high enough and the value 
of Vm in (57) is close enough to unity, then it is not possible to reproduce the 
quantum predictions. In fact, the conditions (9) imply Wm < 1 ( in addition 
to < 1) which give after some algebra the constraint 



V M <\ {V m + V m ) (si + si) + 4^ 



1 - § (si + si) 



, Sj = sm [—) . 



3 ^VlV2 

(59) 

If Tjj « 1 and both V M and V m are close to unity, as is usually the case 
in typical experiments (where V m > 0.96, 77^ < 0.3), the inequality becomes 
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(compare with (28)) 

Vm< + (60) 

I stress that the quantities rjj refer here to the efficiencies of the photon detec- 
tors and not to the overall detection efficiency (which may be substancially 
smaller than rjj due to several losses). Agreement between the model and 
quantum predictions is possible only if eq.(59) (or eq.(60)) is fulfilled. If this 
is not the case, the model predicts inequalities which are violated by quan- 
tum mechanics. In the following I study this case, which allows empirical 
discrimination of the family of local models here proposed versus quantum 
mechanics. 

From the results of section 3 we see that the best choice for the functions 
Pj is provided by eq.(27) (or, what is the same, eq.(58) with /3 - = 1.) "Best" 
is defined in the sense that the disagreement with the quantum predicion is a 
minimum, the disagreement being measured by an appropriate generalization 
of eq.(19). In order to get the best function p we may assume that the 
empirical results to be compared with the LHV model prediction are obtained 
with a fixed 2 ( see ec f(57)). In this case all inequalities of section 4 are valid 
with the effective visibility, V e ff, of the correlation curve substituted for V, 
where we define 

V eff = [V m + (V M - V m ) cos40 2 ] . 

It is easy to see that the stringent inequalities are obtained when V e ff = Vm, 
corresponding to 2 = 0. Thus the inequalities derived in the previous section 
apply, with Vm substituted for V, including the extension to experiments 
using two-channel polarizers and four detectors (a trivial change is required 
if the efficiencies of the detectors are not the same). 

6 Discussion 

Empirical tests of the models defined by eqs.(7) and (8) are possible by means 
of the inequalities (28) and (42), any experiment where one of these inequal- 
ities holds true being compatible with a local hidden variables model of the 
family. The tests are not difficult because visibilities about 97% and detection 
efficiencies of the order of 20% would be enough and these conditions have al- 
ready been achieved in performed experiments. For instance the experiment 
by Kurtsiefer et al.[6] reports V M = 0.982 ±0.001 and V m = 0.970 ±0.001 af- 
ter an appropriate subtraction of accidentals. Any of these values, combined 
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with the efficiency rj = 0.214, leads to a violation of the inequality (28), thus 
making the test possible using the inequalities (42) or (47) . The quantities 
reported in the published paper do not allow a test of the inequalities but the 
experiment clearly shows that the empirical tests may be easily performed 
with present technology. 

The discrimination between the family of local models and quantum me- 
chanics is substantially more difficult than just to test the models. In fact, 
calculating the predictions of quantum mechanics in actual (non-ideal) ex- 
periments is far from trivial. For instance, in the ideal case the quantum 
prediction is given by eq.(17) with both the detection efficiency and the vis- 
ibility of the coincidence curve equal to unity, i. e. V — rj — 1. With these 
values both inequalities (28) and (42) are violated, which excludes all local 
models of the family here studied. This agrees with the well known fact that 
no local hidden variables model is compatible with quantum mechanics for 
all (ideal) experiments. Nevertheless the non-perfect behaviour of detectors 
lowers their efficiency, 77, making the results compatible with local models in 
all experiments performed till now, also a well known fact [2]. 

A restricted family of local models, like the one studied in this paper, may 
be empirically refuted with a moderate value of the detection efficiency, 77, 
combined with a relatively high value of the visibility, V, provided that there 
is rotational symmetry, as assumed in section 3. However in the fashionable 
parametric fluorescence experiments rotational invariance does not hold true 
which makes easier the construction of local models compatible with the 
experiments. In particular it is not obvious that quantum mechanics, with 
all non-idealities taken into account, predicts a detection counting rate of the 
form (57) . If this is not the case the fulfillement of the inequalities (42) or 
(47) would not imply a violation of quantum predictions. 

In summary, the non-ideal behaviour of any actual experiment increases 
substantially the range of parameters where quantum mechanics is compati- 
ble with local hidden variables theories, thus making the empirical discrimi- 
nation rather difficult. It may even be the case that no actual, i. e. non-ideal, 
experiment allows discriminating local hidden variables versus quantum me- 
chanics. 
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